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Abstract 

It is conjectured by Pugh, Lins and de Melo in [7] that the system of equations 

X = y — F{x) 
y = -x 

has at most n limit cycles when the degree of F = 2n + 1 or 2?! + 2. Put M for uniform upper 
bound of the number of limit cycles of all systems of equations of the form 

X = y — {ax"^ + hx'^ + cx^ + dx) 
y = 

In this article, we show that M 7^ 2. In fact, if an example with two limit cycle existed, one could 
give not only an example with n + 2 limit cycles for the first system, but also one could give a 
counterexample to the conjecture A''(2, 3) = 2 [see the conjecture A''(2, 3) = 2 of F. Dumortier and 
C.Li, Quadratic Lienard equation with Quadratic Damping, J. Differential Equation, 139 (1997) 
4159]. We will also pose a question about complete integrability of Hamiltonian systems in R'* 
which naturally arise from planner Lienard equation. Finally, considering the Lienard equation as a 
complex differential equation, we suggest a related problem which is a particular case of conjecture. 
We also observe that the Lienard vector fields have often trivial centralizers among polynomial 
vector fields. 



2000 AMS subject classification: 34C07. 
Keywords: Limit cycles, Lienard equation. 



1 Introduction 



We consider Lienard equation in the form 

\ X = y — F{x) 



(1) 



where F{x) is a polynomial of degree 2n + 1 or 2n + 2. It is conjectured in [7] that this system has at 
most n limit cycles. In particular, it was conjectured that the system 

{X ~ y — {ax'^ + bx^ + cx^ + dx) 
y=-x 

has at most one limit cycle. The phase portrait of (1) is presented in [7]. (1) has a center at the origin 
if and only if F{x) is an even polynomial. The following useful Lemma is proved in [7]: 

Lemma 1 Let F{x) = E{x) + 0{x) where E is an even polynomial and O is an odd polynomial, and 
that 0{x) = has a unique root at x = 0. Then system (1) does not have a closed orbit. 

To prove the Lemma, it was shown that a first integral of the system 

(3) 




that is analytic and defined on — is a monotone function along the solutions of (1). We introduce 
the following conjecture about the above first integral: 



Conjecture 1 Let E{x) be an even polynomial of degree at least then there is no global analytic first 
integral for system (3) on R^. 

The reason for conjecture: There are two candidates for defining a first integral, namely the square of 
the intersection of the solution with negative ?/axis and the other the square of the intersection of the 
solution with .xaxis. The first is well defined on but certainly is not analytic at the origin and the 
second is analytic in the region of closed orbits but it can not be defined on all of R^. In fact, there are 
solutions not intersecting the .xaxis: Consider the region surrounded by y = ^ and j/ = for x 1 
and look at the direction field of 

X = y — F{x) 
y^-x 

on the boundary of this region. We conclude that there is a solution remaining in this area for all t < 
whenever the solution is defined. The analyticity of the second function, intersection with xaxis, is 
obtained by the fact that the Ifrom (y — E{x))dy + xdx is the pull back of (y — g{x))dy + dx under 
TT{x,y) = {x'^,y) with E{x) ~ gix"^)- Now, the differential equation 

x = y- g{x) 

y = -i 

does not have a singular point and the intersection of the orbits with laxis defines an analytic function 
as a first integral, which we call K{x,y). Then K{x^,y) is a first integral for the original system (3). 
Note that putting E{x) = x'^, then the dual from of (3) is the pull back of Ricati equation, whose 
solution can not be determined in terms of elementary functions (This can be seen using Galois Theory, 
see [6]). I thank professor R. Roussarie for hinting me to this pull back. In line of the above conjecture, 
we propose the following question: 
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Question 1 Let an analytic vector field on the plane have a nondegenerate center. As a rule, is it 
possible to define analytic first integrals globally in the domain of periodic solutions? (However, there 
is one locally.) In fact, using Riemann mapping theorem we can assume that the region of the center is 
all ofR^. 

We return to Licnard equation (1). The limit cycles of (1) correspond to fixed points of Poincare return 
map. Let F{x) has odd degree with positive leading coefficient. As a rule, we can not define Poincare 
return map on positive ysemiaxis. For example, if F(x) = a;^"^^ + 2x, in which case the singular point is 
a node, there is no solution starting on positive j/axis and returning again to this axis; see the direction 
field 

x = y- + 2x) 

y = -x 

on the semiline y — x and a; > 0. Then, contrary to what is written in [7] or is pointed out in [10], we 
can define a Poincare return map only in the case that origin is a weak or strong focus or in the case 
of existence of at least one limit cycle. In [10], it is also pointed out that Dulac's problem is trivial for 
(1), that is, for any given F{x), (1) has a finite number of limit cycles. Let F{x) has odd degree, then 
Pq = limy^+oo P{y) exists and even we have \iniy^+oo{P{y) — Po)y'' = for all i. Then P has a finite 
number of fixed points. (In this work we do not use the strong approach limy_,+oo y^{P{y) ^ Pq) ~ 
and thus we do not prove it.) For the case that F has even degree, the above is not trivial and one can 
deduce it from the results of this paper. In fact we must consider the case that limj,_,+oo P{y) = +oo. 
Equivalently, we have a loop at Poincare sphere based at infinity. As a simple consequence of the above 
Lemma, we note that the system 

X ~ y — x^ 
y ~ ^{o- — x) 

does not has limit cycles for a 7^ 0, because putting x := x + a, y :~ y — we will obtain 

X = y — [x"^ + 2ax) 

y = -<^x. 

Using the Lemma, this system does not have a limit cycle. In figures on page 478 and 479 of [4], it 
appears that the existence of limit cycles is claimed. 

2 Main Results 

Theorem 1 For any a, b, c, there exists a unique d ~ (a, b, c) such that the system (2) has a homoclinic 
loop in Poincare sphere. This loop is stable if and only if the singular point is unstable. 

Corollary 1 The maximum number of limit cycles of (2) can not be exactly two. 

Proof: (Proof of the theorem) For fixed a, b, c, a > 0, 6 > 0, put U{d) and S{d) for intersection of 
unstable and stable manifolds corresponding to the topological saddle (0, 1, 0) on the equator of Poincare 
sphere. In fact U{d) and S{d) are similar to P- and P+ on page 339 of [7], figure 3. U{d) and S{d) are 
continuous and monotone functions (as will be proved). We also prove U{d) = S{d) has a unique root. 
First note that for d > there is no limit cycle or homoclinic loop based at infinity, using the lemma. 
Then we assume d < 0, note that U{Q) > S{Q), otherwise from the stability of the origin in (2) for 6 > 
and d = 0, we would have a limit cycle for (2) which is a contradiction with the lemma. On the other 
hand, U{d) ^ S{d) for d <^ —1 because the minimum value of F{x) = ax'^ + bx^ + cx"^ + dx goes to —00 
as d — > —00 and S is decreasing. Now by continuity and monotonicity of U and S we have a unique root 
do = 'ip{a,b,c,) such that U{dQ) — S{do). We must prove that U and S are continuous and monotone. 
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It can be directly observed, even without using any classical theorem, that U, S are continuous. We 
prove that S is decreasing, similarly U is increasing. Let 7^ be a solution of {2)d that is asymptotic to 



the graph of F{x) = ax + bx 



■ cx 



- da:, in a; < in fact, 7^ is the stable manifold for the topological 



saddle (0, 1, 0) that intersects yaxis in S{d). We prove that for d' < d < 0, S{d') > S{d). is a curve 
without contact for (2)^' and the direction field of (2)^^ on 7^ is toward "left" and the direction of (2)^/ 
on the semifine y > 0- x = 0, is toward "right". Then there is a unique orbit of {2)di that remains in 
the region surrounded by 7^ and the semifine x ^ Q, y > 0. This orbit is a stable manifold for (0, 1, 0) 
of system (2)^', say; jd' certainly 7^ can not intersect 7^ so 7^' will intersect negative yaxis in S{d') 
above S{d). Therefore S is decreasing. For the proof of the theorem it remains to prove that the loop 
is attractive, still assuming a > 0, 5 > 0, d < 0. Before proving that the homoclinic loop is attractive 
we point out that, at first glance, this loop has a degenerate vertex at (0, 1,0), i.e. the linear part of 

Q Q . But using weighted compactification as explained in [3], one fias 

an elementary polycyclc with two vertices, one vertex is a " From" vertex and another is a "To" . Thus 
we do not have an "unbalanced polycycle" and therefore behavior of solution near the polycycle can 
not be easily determined. Thus, we use the following direct computation, (for definition of balanced 
and unbalanced polycycles, sec [5] page 21.) Let Pq = U{do) = S{do) < 0; We have a Poincare return 
map P defined on negative vertical section (Pq, 0), we will prove that liniy^p,-, P'(y) = 0, then P{y) — y 
is negative for y near Pq. Note that considering the orbit of points near Pq in the vertical section 
(Po,0) is equivalent to considering the orbits of points (0,y) with y ^ 1. Let 7 be the orbit of (2) 
corresponding to a homoclinic loop whose existence is proved above. 7 is asymptotic to the graph of 



vector field at the vertex is 



Fix) 



bx"^ 



dx. 




Figure 1: F{x) = ax^ + bx^ + cx^ + dx 

The orbit starting from (0,y) intersects the graph of F at points with ^coordinates yi and j/2, 
when traced in positive and negative time direction, respectively. Recall that we want to prove 
liiTiy^Pg P'{y) = 0; we use the following three facts: 

I) There is a constant K such that jy — < and |j/ — 2/2! < Ky/y, where K depends only on 
a, 6, c, d. 

II) 7 is asymptotic to the graph of F. 

III) Let A{y) and B{y) two right inverse of F{x) = ax^ + bx^ ^-cx"^ + dx, that is F{A{y)) = F{B(y)) = y 
for y ^ 1. Then limy^+oc{A{y) + B{y)) = (III) is a simple exercise (II) is pointed out in [7]: So we 
only prove the first. It suffices to prove (I), with the same notation yi and j/2 for the intersection of the 
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orbit with the graph of F{x) + 1 (in place of F{x)).\y — yi| = -l^^^Ax, where Ax is the ^coordinate of 
the intersection of the orbit starting from (0, y) with the graph of F{x) + 1, then Ax < K^/y for some 
K, since the degree of F is 4 and = y^Sfl^a) ' '^'^'5'"*^ (^O: Uo) is a point of the trajectory starting 
at (0,y) and < a; < Aa; (using mean value theorem). Then ^^^r < K^/y and (I) is proved. Now we 
compute liniy^Pg P'{y) (or cquivalcntly for y ^ 1 as stated above): P'{y) = 'pfyj^^Pilo^^'' ''div"dt] 
where by "div" we mean the divergence of vector field (2). i.e. —F'{x{t)), and T{y) is the time of first 
return of the solution starting {Q,y) to negative vertical section (for the standard formula of P'{y) see 
[I])- 

' , , /■ X , f F'(x) f -F'(x)dr 

divdt = - / F'{x)dt ^ ' ^ ' ' ^ ' 



y~F{x) 

We use y for the intersection of the solution starting at (0, y) with the positive j/axis in positive time. 
We show that /q^^^^ div = J^^^^ —F'{x{t))dt goes to —oo for y near Pq. We divide /q^^^^ divdt into 
three parts /i, I2, Is'- Ii for the part of the integral that the orbit jy, the solution starting at {0,y) for 
y G (pojO) lice in \x\ < c, where c is a nonzero constant. I2 corresponds to the part of above the 
horizontal line y — y — K's/y where K is the same constant which is given in (I), and /a is the remaining 
part of / = /o'"^^^ fact, in Iz wc compute the integral of the divergence of (2) along the part of 

7y that lies below the horizontal line y — y ~ Ky/y and outside of \x\ < ci where in this part jy lies 
between the graph of F{x) and the orbit 7 which corresponds to the homoclinic loop. Note that as 
y +00 and consequently y Pq, I2 and I3 are very large in absolute value but Ii remains bounded 
because / F'{x) = J ^-rf^- On the other hand we will see that I2 +I3 goes to — 00 as y ^ Pqi thus we 
may ignore the term Ii. We could choose K in (I) so that |-^— < KVv where A{y) is the positive 

inverse of F{x) for large y. There is a constant Ki such that I/2I < Kiy. Now we show to —00 as 
y ^ Pq. Recall that y is as in the figure. Applying III, we realize that I3 is the integral of some function 
which is big in absolute value as y Pq (or y +00). Generally speaking, let (7(7) be a function such 
that limj,^+oo 3(7) = +00. Put G(7) = Q g{s)ds, then limy^_|_oo = +00. So, it suffices to prove 
that I3 is the integral of some function with respect to y, where this function goes to —00 as y —f +00. 

,,^f^.l,. (4) 



F'(x) 

— 402;^ + 3&.T + 2c + dx. 

X 

We consider two parts of 7^, one in x > c and another in a; < — ci since below the horizontal line 
y = y — Ky/y the orbit 7j, lies between 7 and graph of F{x). Note that 7 is asymptotic to the graph 
of F{x). Now, applying III, we will obtain — —00: for large values of yi the term in (1) can 
be omitted, then we must compute — hnij^+oo 4a(A^(y) — B^{y)) + 3b{A{y) — B{y)), where A{y) and 
B{y) are inverses of F{x) such that F(A(y)) = F{B{y)) = y and A{y) > 0, B{y) < 0. We look at 
-limy+oo(^(y) - B{y))[4:a{A{y) - B[y)) + 36]. Certainly {A{y) - B{y)) +00. So, the above limit 
goes to —00 and +00 as y ^ Pq- This completes the proof of the theorem. □ 

Proof: (Proof of the corollary) From the proof of the theorem, we conclude that (2) can have an even 
number of limit cycles for d < do ~ 4' {a, b, c) and can have odd number of limit cycles for do < d < 0. 
Now let {2)d have exactly two limit cycles, therefore d < do = ip{a,b,c). If d < do, then (2)^^ has at 
least two limit cycles, counting multiplicity. It is because any closed orbit of (2)^ is a curve without 
contact for {2)cto and the direction fields of (2)^^^ on this closed curve are toward the interior. But {2)cig 
can not have an odd number of limit cycles. For instance let none of the two limit cycles of (2)^;^ be 
semistable. Then by a small perturbation d = dg — e, (2); these two limit cycles do not die. On the other 
hand, we would obtain another limit cycle near the loop 7. Because for d = do — e, 7 is a curve without 
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contact for (2)^ and the direction of (2)^ on 7 is toward the interior. Note that 7 was an attractive 
loop, therefore if e is very smah, using PoincarcBendixson Theorem we will obtain a third limit cycle 
near the loop. Now any semistable limit cycle can be replaced by two limit cycles by an appropriate 
perturbation d = doe depending on which side of the limit cycle is stable. For instant, let (2)^;, have a 
semistable limit cycle in interior of the loop then d = do — e gives us two limit cycle near the semistable 
ones and simultaneously one limit cycle near the loop. However the corollary is proved, we point out 
that the existence of 3 limit cycles for (2) easily implies the existence of 3 limit cycles for 



for small e. Now by a linear change of coordinates we put e = 1, so we would obtain counterexample to 
the conjecture that the latter system for e = 1 has at most 2 limit cycles. See conjecture A^(2, 3) = 2 in 



Remark 1. For the proof of the existence of the loop 7 and also the proof of the corollary, in fact we 
used the rotational property of the parameter d in (2), namely, any solution of (2)^ is a curve without 
contact for (2)^'. if d' ^ d. In particular, periodic solutions of {2)d are closed curves without contact for 
{2)d'. The simple but useful phenomenon of "rotated vector field theory" introduced by Duff, is some 
times used erroneously. See, for example , the investigation of 



in [8]. It is claimed that there is no limit cycle for d < do ~ -(/'(a, 6, c) < 0, where 6 > 0. In fact, the 
following situation that could occur, is not considered in [8]: 

As we assumed above, let a > 0, b > 0, for small c? < we have exactly one (small) Hopf bifurcat- 
ing limit cycle. It is possible that this limit cycle, before arriving to a loop situation, dies out in a 
semistable limit cycle. Put X = (i|(5)(;hasexactlyonelimitcycle, we do not necessarily have do = infX, 
where do = ^'(Oj^jO), correspond to loop situation. Put "i" for the above infimum. It could be i > do 
and (5)i possesses a semistable limit cycle. It is also possible that when the outermost limit cycle is 
dying out in the loop, the two innermost limit cycles have not arrived to each other yet. In fact [8] 
suggests an affirmative answer to the conjecture for system (2). 

Remark 2. The homoclinic loop 7, as an orbit on the plane, not on the Poincare sphere, divides 
the plane into two parts: its interior, where all solutions are complete, and its exterior, where all solu- 
tions have a finite interval of definition. Interior orbits are complete because 7 is a complete orbit by 
virtue of J dt ~ J ^J^^^^^ and 7 being asymptotic to the graph of F{x). The exterior points of 7 are 
not complete orbits because they tend to hyperbolic sink and the source on the equator of the Poincare 
sphere (See [2]). Now, a trivial observation is that Licnard equation (1) can not have an isochronous 
center i.e. a center with a fixed period for all closed orbits surrounding it. 

Remark 3. As we saw in the proof of the theorem and the corollary, the inequality U{d) < S{d) 
or the reverse, determines oddness or evenness of the number of limit cycle. In this direction we point 
out that: Let F{x) be an even degree polynomial with positive leading coefficient, and U{e) and 5(e) 
be similar to U{d) and S{d) above for the system 



Then lim^^o U{e) ~ A/+ and lime_,o S{e) = M where M and are minimum values of F{x) on 
(—00, 0] and [0, +00) resp. The proof is identical to the proof of existence of orbits passing through 




[3]. 



□ 
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U{d) and S{d) asymptotic to the graph of F{x), as in [7]. 

Remark 4. Note that when the degree of F{x) in (1) is odd, the behavior of infinity is determined 
only by the sign of the leading term of F{x). Then, giving an example of (2) with 3 limit cycles would 
give, inductively, n + 2 limit cycles for (1), for all n. 

Remark 5. Considering "flow" version of the problem of "centralizer of diffeomorphisms" described 
in [10] one can easily observe the following partial result. Let L be the Licnard vector field similar 
to (1) with at least one closed orbit and X be a polynomial vector field such that [i,^] = 0, then 
X ~ cL where c is a constant real number. In general, let two vector fields have commuting flows and 
7 be a closed orbit for one of them which does not lie in an isochronous band of closed orbits. Then 
7 must be invariant by another vector field and if both vector fields are polynomials. Then either 7 
is an algebraic curve or two vector fields are constant multiple of each other. But Lienard systems do 
not have algebraic solutions, [9]. More generally by the following proposition we have "Non existence 
of algebraic solution implies triviality of centralizer" . 

Proposition 1 Let M be the set of all polynomial vector fields on C^. Then X E M has trivial 
centralizer if X dose not have an algebraic solution. 

Proof: Let [X,Y] = 0. Then X.Det{X,Y) = {DivX).{Det{X,Y)), so ''Det{X,Y) ^ 0" defines an 
algebraic curve invariant under X (By Det{X, Y)) we mean the determinant of a 22 matrix whose 
columns are components of X and Y). □ 

Remark 6. The following could be a (real) generalization of formula used in proof of proposition: 
Let (M, uj) be a 2n dimensional symplectic manifold and X, Y two vector fields with the condition 
[X, y] = 0, then X.lo{X,Y) = n{DivX)uj{X,Y). This formula is trivial for usual symplectic structure 
of . Further there is a local chart around each point of a two dimensional symplectic manifold that 
to can be represented in the trivial form. Thus the formula is proved for arbitrary two dimensional 
symplectic manifold. Now, in general case we have a two dimensional symplectic submanifold N of A/, 
that X and Y are tangent to N. (For points m G M that uj{X{m), Y{m)) 7^ 0, using Frobenius theorem.) 
From other hand DivXM — nDivXpf. The investigation of points m that uj(X(m), Y(m)) = is trivial. 
Then the proof is completed. 

Question 2 Dirac introduced the following embedding of planner vector fields into Hamiltonian system 
inW^:H = zP{x, y) + wQ{x, y). Now we consider the Hamiltonian H — z{y — F(x)) — wx. When is 
this Hamiltonian completely integrable? 

By completelyintegrable Hamiltonian, we mean that there is a first integral for the system 

y - Fix) 

— X 

w — F'{x)z 

—z 

independent of H . 

The particular case F(x) = x'^ for which (1) has a center with global first integral (j){x,y) = {y — 
x^ + ^)e~'^y , suggests that when (1) does not have a center, the above Hamiltonian is not completely 
integrable. In fact, the function (j), as above, is a first integral, independent of H for above four 
dimensional system. Then, when (1) is not integrable, one expects that the corresponding Hamiltonian 
is not completely integrable. However, surprisingly, putting F{x) = kx, we have another first integral 
yw + xz. 
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Question 3 Considering (1) as a vector field on C^, and in line of conjecture in [7] one can think of 
the validity of the following two statement. 

I) There are at most n different leaves containing real limit cycles. 

II) There are at most n real limit cycles lying on the same leaf. By "leaf" we mean a leaf of the foliation 
corresponding to the equation (1) on — 0. 

Acknowledgement. The author thanks Professor S. Shahshahani for fruitful conversations. He 
also acknowledges the financial support of the Institute for Studies in Theoretical Physics and Mathe- 
matics. 
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